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Abstract 

A search for the lightest neutral CP-even and neutral CP-odd Higgs bosons 
of the Minimal Supersymmetric Standard Model is performed using 216.6 plT 1 of 
data collected with the L3 detector at LEP at centre-of-mass energies between 203 
and 209 GeV. No indication of a signal is found. Including our results from lower 
centre-of-mass energies, lower limits on the Higgs boson masses are set as a function 
of tan/3 for several scenarios. For tan/5 greater than 0.7 they are > 84.5 GeV 
and rriA > 86.3 GeV at 95% confidence level. 
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1 Introduction 



In the Minimal Supersymmetric Standard Model (MSSM) pQ two doublets of complex scalar 
fields are required to generate the masses of gauge bosons and fermions. The neutral Higgs 
sector of the MSSM comprises three physical states: two CP-even Higgs bosons, the lighter of 
which is denoted as h and the heavier as H, and a neutral CP-odd boson, A. 

The two most important production mechanisms of the light neutral Higgs boson in e + e~ 
collisions are: 

e+e" -> hZ, (1) 



e + e" -> hA, (2) 

with tree level cross sections that are related to the Standard Model Higgs-strahlung cross 
section, cr^z, as 0: 

a hZ = sin 2 (/?-a)4^, (3) 



a hA = cos 2 ([3- a)~\a™, (4) 

where tan/? is the ratio of the vacuum expectation values of the two Higgs doublets, a is 
the mixing angle in the CP-even Higgs boson sector and A is the p-wave suppression factor 
depending on the Higgs boson masses, and m^, and the centre-of-mass energy y/s. For some 
choices of the MSSM parameters, the H boson can also be produced via the Higgs-strahlung 
process 

e+e" -»■ HZ, (5) 

with a cross section suppressed by the factor cos 2 (/3 — a) relative to cr™ . 

For most of the MSSM parameter space considered, the neutral Higgs bosons are predicted 
to decay dominantly into bb and t + t~. However, in certain parameter regions, other decays 
like h — > AA and A —>■ cc become important. 

The search for the neutral Higgs bosons is performed in the framework of the constrained 
MSSM with seven free parameters. These are the universal sfermion mass parameter, Msusy, 
the common Higgs-squark trilinear coupling, A, the supersymmetric Higgs mass parameter, /i, 
the SU(2) gaugino mass parameter, M 2 , the gluino mass parameter, M 3 , m A and tan/?. The 
mass of the top quark is taken to be 174.3 GeV j^j. 

Previous searches for the neutral Higgs bosons were reported by L3 [HE] and other experi- 
ments jH] . In this Letter, we present the results of the search for the h and A bosons using data 
collected with the L3 detector [Z| in the year 2000. In comparison to previous analyses [I] the 
performance in the four-jet channel is improved by using a new likelihood based analysis. In 
addition, the six-jet final state resulting from the hZ — > AAqq topology is investigated. 



2 Benchmark Scenarios 

Due to the large number of parameters remaining in the constrained MSSM, we focus on 
three specific parameter settings as suggested in Reference |S1 These benchmark settings are 



2 



denoted as the "mh— max", "no mixing" and "large— /i" scenarios, with the corresponding 
MSSM parameter values detailed in Tabled The first two scenarios take into account radiative 
corrections to computed within a two-loop diagrammatic approach |Hj and differ only by the 
value of X t = A — ficotft, which governs the mixing in the scalar top sector. The "mh— max" 
scenario is designed to extend the search to the maximal theoretical bound on for any 
value of tan/3 and leads to rather conservative constraints on tan/3. The "no mixing" scenario 
corresponds to vanishing mixing in the scalar top sector and is more favorable to LEP searches. 
After fixing the parameters Msusy, M2, /i, X t and M3, a scan over the remaining parameters 
tan/? and is performed in the range 0.4 < tan/3 < 30 and 10 GeV < < 1 TeV. The 
widths of the h and A bosons are assumed to be smaller than the experimental mass resolutions. 
As this holds in these two scenarios only for tan/3 < 30, higher values of tan/3 are not considered. 

In the "large— /x" scenario the upper theoretical bound on m-h is slightly less than 108 GeV, 
thus in the LEP reach. However, for some choices of tan/3 and m-A, the Higgs boson pair 
production (j2J) is kinematically inaccessible and the Higgs-strahlung process is strongly 
suppressed due to small values of sin 2 (/? — a). Sensitivity in these regions can be recovered 
by exploiting the H boson production via the high cross section Higgs-strahlung process (0). 
The scan is performed over from 10 GeV to 400 GeV and over tan/3 between 0.7 and 50. 
Within this scenario, the assumption that the widths of the Higgs bosons are small compared 
to experimental mass resolutions is valid for tan/3 up to 50. For the interpretation of the data 
within this "large— /z" scenario, the improved one-loop renormalisation group calculations JU] 
are used. 

3 Data and Monte Carlo Samples 

The data recorded at centre-of-mass energies between 203 and 209 GeV are grouped into three 
data sets with effective centre-of-mass energies of 204.3, 206.1 and 208.0 GeV, corresponding 
to integrated luminosities of 26.0, 181.9 and 8.7 pb~ . The results obtained from this data 
are combined with the results from integrated luminosities of 233.2 pb _1 at 192 GeV < y/s < 
202 GeV 4J and 176.4 pb" 1 at y/a = 189 GeV 0. 

The cross sections of processes ([Q, (0) and (j3j) and the decay branching fractions of h, 
H and A are calculated using the HZHA generator |11| . For efficiency studies, Monte Carlo 
samples of 2000 Higgs events are generated for each mass hypothesis in each search channel 
using PYTHIA jT2j and HZHA. For the hA samples, mh and rriA range from 50 to 105 GeV in 
steps of 5 GeV. For the hZ samples, m h is chosen in steps of 5 GeV from 60 to 100 GeV and 
in steps of 1 GeV from 100 to 120 GeV. For background studies, the following Monte Carlo 
programs are used: KK2f [T3] (e + e~ — > qq(7)), PYTHIA (e + e~ — > ZZ and e + e~ — > Ze + e~), 
KORALW HU (e + e~ -> W+W") and KORALZ (TH] (e+e" -> r+r"). Hadron production in 
two-photon interactions is simulated with PYTHIA and PHOJET jTEj. EXCALIBUR [Hj is 
used for other four fermion final states. The number of simulated events for the most important 
background channels is more than 100 times the number of expected events. 

The L3 detector response is simulated using the GEANT program ^H]; which models the 
effects of energy loss, multiple scattering and showering in the detector. The GHEISHA pro- 
gram [THj is used to simulate hadronic interactions in the detector. Time dependent detector 
inefficiencies, monitored during data taking, are also taken into account. 
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4 Analysis Procedures 



For the hA production, the decay modes considered are: hA — > bbbb, hA — > bbr + T~ and 
hA — > r + r _ bb. For the hZ or HZ production, the event topologies considered are: qqq'q', 
qqvv, qq£ + £~(£ = e, //, r) and r + r~qq. To cover the MSSM parameter regions where the decay 
h — >• AA becomes important, the hZ — > AAqq — > qqq'q' q"q" channel is studied as well. Searches 
in channels with decays of the h boson into quarks are optimised for the dominant h — > bb de- 
cay channel. The analyses of the qquu and qq£ + £~(£ = e,/i) final states are the same as those 
used in the Standard Model Higgs search [2U] . 



4.1 The hZ — > bbqq and hA — > bbbb Analyses 

The signature of the hZ — ► bbqq final state is four high-multiplicity hadronic jets and the 
presence of b hadrons, in particular in the jets expected to stem from the Higgs boson. The 
invariant mass of the jets supposed to originate from the Z boson must be compatible, within 
mass resolution, with its mass, mz- The hA — > bbbb search topology is characterised by four 
high multiplicity hadronic jets originating from b-quarks. The main backgrounds arise from 
qq(7) final states and hadronic decays of W-pairs and Z-pairs. 

Initially, a common preselection for both the hZ — > bbqq and hA — > bbbb channels is ap- 
plied, followed by a kinematic classification of events into hZ or hA analysis branches. Selection 
criteria optimised for each branch are used and a final discriminant specific to each branch is 
constructed. 

The preselection criteria used are described in Reference EHI Events passing the preselec- 
tion are forced into a four-jet topology using the Durham algorithm [21] and a kinematic fit 
imposing energy and momentum conservation (4C fit) is performed. Each event is tested for 
its compatibility with the hZ and hA production hypotheses, exploiting dijet invariant masses. 
There are three possibilities for pairing jets in a four-jet event. For each pairing, \ 2 values are 
calculated for the hypotheses of hZ and hA production: 

2 (Si - m h - m z f (Aj - \m h - m z |) 2 

Xhz — ~2 1 2 ' W 

a S °A 



2 (Si -m h -m A ) , (Ai-|m h -m A |) 

XhA = ~2 + ~2 • (7) 

In these expressions Sj and Aj are the dijet mass sum and dijet mass difference of the i-th pair- 
ing, and as and a a are the dijet mass sum and dijet mass difference resolution functions. They 
are found to be almost independent of the dijet masses for production well above threshold and 
are estimated from Monte Carlo to be = 4 GeV and oa = 18 GeV. Close to the kinematic 
threshold they strongly depend on the dijet mass sum. For each hypothesis, the jet pairing 
with the smallest \ 2 is chosen and the probability P(x 2 ) is calculated. 

Events are assigned either to the hZ or hA analysis branch by means of a binned likeli- 
hood [213, L£ z = 1 — L£ A , constructed from the following variables: the x 2 probabilities, P(Xhz) 
and P(XhA); ^ ne Higgs boson production angle, | cos9|, the number of charged tracks, N trk , the 
global event b-tag, B tag , and the maximum triple jet boost, 7 tr ipie |2IJ|- Events are assigned to 
the hZ branch if L£ z > 0.5, or to the hA branch otherwise. 

The selection in both branches proceeds in the same way. Events are rejected if P(x 2 ) < !%• 
High b-tag events are accepted and a selection likelihood is then constructed to separate the 
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signal from the qq(7) final states and hadronic decays of W-pairs and Z-pairs, using N tr k, 
Btag, Ttripiej I cos6| and logY 34 , where Y 34 is the jet resolution parameter for which the event 
topology changes from three to four jets. A final discriminant is constructed for events passing 
an optimised selection likelihood cut. The optimisation is based on the analysis performance 
at (mAjfRh) values close to the expected sensitivity of the L3 combined search. The final 
discriminants are built from individual b-tag variables of the four jets and P(x 2 )- 111 the hZ 
branch, an event category variable is also used. This variable is constructed from rankings of 
b-tag variables of the two jets assigned to the Higgs boson, as described in Reference EDI The 
final discriminant in the hA branch also includes the corresponding selection likelihood. 

Table[2]reports the number of data, expected background and expected signal events selected 
at different stages of the analysis for two representative Higgs boson mass hypotheses in the 
"mh— max" scenario, {m^rn^) = (90,90) GeV at tan/3 = 25 and (m^m^) = (165,110) GeV at 
tan/5 = 3. For these hypotheses the distribution of the final discriminant in terms of a signal- 
to-background ratio is presented in Figure ^ The data are compatible with the background 
expectation. 

4.2 The hZ — > bbr + r _ , r + r _ qq and hA — > bbr + r _ Analyses 

The signatures of hZ — > bbr + T~, r + r _ qq and hA — > bbr + T~ final states are a pair of taus 
accompanied by two hadronic jets. For each of the channels hA and hZ an analysis is optimised 
based either on the tau identification or on the event topology by requiring four jets with two 
of them being narrow and of low multiplicity. The main background comes from W-pair decays 
containing taus. The analysis is similar to the one used in previous searches jlj. The selection 
is optimised for lower Higgs boson masses by applying looser cuts on the opening angles of the 
jets and tau pairs compared to the Standard Model Higgs search [20]. The invariant mass of 
the tau pair, m T + T -, and of the hadronic jets, m q q, must be between 25 GeV and 125 GeV. The 
ratio between the total energy deposited in the detector, E vis , and yfs must be less than 0.9 
and the polar angle of the missing momentum, O m i SS , must satisfy | cos0 m i SS | < 0.9. 

Four possible final states are considered: hZ — > bbr + r~, hZ — > r + r~qq, hA — > bbr + r~ and 
hA — ■> r + r~bb. Each event is uniquely assigned to one channel using mass and b-tag informa- 
tion, and final discriminants specific to each channel are constructed. For the hZ — > bbr + T~, 
hA — > bbr + T~ and hA — ► r + r~bb final states, the final discriminant is constructed from m qq , 
m T + T - and the b-tag variables of the two hadronic jets. For the hZ — > r + T~qq final state, 
m T + T - is used as the final discriminant. 

Table El reports the number of data, expected background and expected signal events for the 
same Higgs boson masses chosen in the previous section. For these hypotheses the distribution 
of the final discriminant in terms of the signal-to-background ratio is shown in Figure 121 Good 
agreement between the data and the expected background is found. 

4.3 The hZ — ► AAqq Channel 

To improve the search sensitivity in the region of low tan/3 and low where the h — > AA decay 
becomes dominant and the A — > cc decay replaces A — > bb, a dedicated analysis is devised 
and performed on the data collected at y/s = 189 — 209 GeV. This analysis aims to select 
hZ — > AAqq — ► qqq'q / q"q" final states and is derived from the analysis used in the four-jet 
channel. At the first stage, the same preselection as in the four-jet channel is applied with 

^Both of the decay modes (h— ► bb, A— > t+t - ) and (h— > t+t - , A— > bb) are considered. 
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an additional cut on the event thrust, T < 0.9. In the next step, a signal likelihood LAAqq is 
built to distinguish the hZ — ► AAqq signal from the qq(7) final states and hadronic decays of 
W-pairs and Z-pairs. This likelihood is constructed from the variables: N tr k, Ttripiej logY34, the 
event sphericity, the absolute value of the cosine of the polar production angle, assuming the 
production of a pair of gauge bosons |20] and logY 56 , where Y 56 is the jet resolution parameter 
for which the event topology changes from five to six jets. Among these variables N tr k and logYs6 
have the most discriminating power between the hZ — > AAqq signal and four-fermion and two- 
fermion backgrounds. The likelihood LAAqq is used as the final discriminant. No evidence for 
the hZ — > AAqq signal is found in data. As an example, Figure |3] shows the distributions of 
logY 56 and LAAqq for data, the expected background and the signal corresponding to the Higgs 
boson mass hypothesis (m^m^) = (30,70) GeV. 

5 Results 

The analyses presented in this Letter are combined with the hZ — > bbvv and hZ — > bb£ + £~ 
(£ = e, n) analyses used in the Standard Model Higgs searches [20] . The results of previous 
searches at lower y/s are also included. The final discriminant distributions obtained in 
each search channel at each centre-of-mass energy are used to evaluate the presence of a signal in 
the data. No evidence for a signal is found and the search results are interpreted in terms of an 
exclusion of MSSM parameter regions. The statistical procedure adopted for the interpretation 
of the data and the definition of the confidence level CL S are described in Reference 1221 The 
analysis performance is quantified with the expected median confidence level, CL me d, which is 
obtained from CL S by replacing the observed value of the test-statistic — 21nQ by its background 
median value. 

Systematic and statistical uncertainties on the signal and on the background are incorpo- 
rated in the confidence level calculations as described in Reference The main sources of 
systematic uncertainties are the detector resolution, the selection procedure, theoretical uncer- 
tainties and Monte Carlo statistics. The overall systematic uncertainty is estimated to range 
from 3% to 6% on the expected signal, and from 7% to 15% on the background, depending on 
the search channel. 

5.1 Limits in the "rah— max" and "No mixing" Scenarios 

Figure |U shows the area of the (tan/3, m^) and (tan/5,mA) planes excluded at 95% confidence 
level for the "mh- max" and "no mixing" scenarios. In the "mh— max" scenario, lower limits 
on the masses of the h and A bosons are set at 95% confidence level as: 

m A > 86.5 GeV, m h > 86.0 GeV, 

for every tan/? value considered. The expected values in the absence of a signal are m\ > 
88.6 GeV and > 88.4 GeV. For 0.55 < tan/5 < 2.2 the A boson is excluded up to a mass of 
1 TeV thus allowing to rule out this tan/5 range. 

In the "no mixing" scenario, the combined results establish lower mass bounds at 95% 
confidence level of: 

m A > 86.3 GeV, m h > 85.5 GeV. 

The expected limits are m A > 88.6 GeV and > 88.5 GeV. Here the tan/5 range between 0.4 
and 5.4 is excluded at 95% confidence level, also for m A up to 1 TeV. A downward fluctuation 
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of about la compared to the background expectation is observed in the data at {mxiTrih) ~ 
(90,90) GeV. There is a deficit of candidates in the hA — > bbbb channel at -^/s = 203 — 209 GeV 
resulting in an observation of a CL S value smaller than 5% although no exclusion of this region 
at 95% confidence level is expected. This effect explains the irregularity in the exclusion plots 
at high tan/? and between 86 GeV and 91 GeV. An excess in the hZ — > bb£ + £~(£ = e, ^) 
channel at ~ 90 GeV and in the hZ — > bbqq channel in the range between 90 GeV and 
95 GeV, from the a/s = 192 — 202 GeV data, result in a sizable reduction of the excluded range 
of tan/3 for the mass range 90 GeV < < 100 GeV in both scenarios. The area at tan/3 < 0.8 
and mA < 40 GeV which previously was not excluded in the "no mixing" scenario [I] is now 
excluded using the results of the hZ — > AAqq analysi^- 



5.2 Limits in the "Large— ^" Scenario 

In the "large— /i" scenario there are regions of the (tan/?,mA) plane where the hA process is 
inaccessible and the hZ process is suppressed by a small value of sin 2 (/3 — a). However, the 
heavy CP even Higgs boson H is expected to be produced there via the Higgs-strahlung process. 
Hence, the loss of the sensitivity for the h boson can be compensated by reinterpreting the hZ 
analyses in the context of the HZ search. This is illustrated in Figure which presents the 
(1 — CL mc d) confidence level calculated as a function of rriA at tan/3 = 15 in the context of 
the searches for the h and H bosons. Searches for the h boson alone lack sensitivity in the 
range 89 GeV < ttia ;$ 108 GeV. The inclusion of the HZ search results extends the region 
of sensitivity leaving only the range 89 GeV < < 97 GeV unexcluded. However, since all 
the analyses, except for the hZ — > r + r~qq channel, are optimised for the Higgs boson decays 
into bb, they do not provide sufficient sensitivity to the parameter regions where the effective 
couplings Hbb and hbb are suppressed [24] . The observed exclusion is presented in Figure Eb- 
Exclusion plots in the (tan/?,mh) and (tan/?,mA) planes for the "large— /z" scenario are 
presented in Figures Eb and Eli, respectively. Limits on Higgs boson masses are derived as: 

m h > 84.5 GeV, m A > 86.5 GeV. 

The expected values are 87.2 GeV and 89.2 GeV, respectively. Furthermore, the range 0.7 < 
tan/? < 6.7 is excluded, for values of up to 400 GeV. The allowed area between tan/? =15 
and tan/3 = 50 corresponds to reduced Hbb or hbb couplings. The unexcluded region at 
~ 88 GeV and tan/3 ~ 15 is caused by a slight upward fluctuation in the data, coming mainly 
from hA — > bbbb candidates selected in the data at > 203 GeV. The allowed vertical narrow 
band at = 107 — 110 GeV and tan/? > 10 represents the region where the hA production 
is kinematically inaccessible and cos 2 (/3 — a) ~ sin 2 (/3 — a) ~ 0.5 so that both the hZ and 
HZ production cross sections are reduced by a factor of 2 with respect to u^iz- Although 
the L3 combined search has a sensitivity for exclusion of this critical region, as can be seen 
from Figure Eb ; the expected median confidence level is only slightly lower than 5% and an 
insignificant upward fluctuation observed in the data pushes the observed confidence level, CL S , 
above 5%, thus not allowing to exclude this region at 95% confidence level. Finally, the allowed 
area at tan/3 between 6.7 and 10 and between 90 and 100 GeV arises due to the excesses 
already discussed. 

In conclusion, no evidence for neutral Higgs bosons of the MSSM is found and large regions 
of its parameter space are excluded. 

2 ^The hZ —> AAqq analysis is used instead of the four-jet one whenever it provides better sensitivity, i.e. 
gives smaller values of CL mec i . 



7 



References 

[1] H. P. Nilles, Phys. Rep. 110 (1984) 1;H. E. Haber and G. L. Kane, Phys. Rep. 117 (1985) 
75;R. Barbieri, Riv. Nuovo Cim. 11 n°4 (1988) 1 

[2] J.F. Gunion et al, in The Higgs Hunter's Guide, (Addison Wesley, 1990), p. 191 

[3] D.E. Groom et al, Eur. Phys. J. C 15 (2000) 1 

[4] L3 Collaboration, M. Acciarri et al, Phys. Lett. B 503 (2001) 21 

[5] L3 Collaboration, M. Acciarri et al, Phys. Lett. B 471 (1999) 321 

[6] OPAL Collaboration, G. Abbiendi et al, Eur. Phys. J. C 12 (2000) 567;ALEPH Col- 
laboration, A. Heister et al, Phys. Lett. B 526 (2002) 191;DELPHI Collaboration, J. 
Abdallah et al, Eur. Phys. J. C 23 (2002) 409 

[7] L3 Collab., B. Adeva et al, Nucl. Instr. Meth. A 289 (1990) 35; 
J. A. Bakken et al, Nucl. Instr. Meth. A 275 (1989) 81; 
O. Adriani et al, Nucl. Instr. Meth. A 302 (1991) 53; 
B. Adeva et al, Nucl. Instr. Meth. A 323 (1992) 109; 
K. Deiters et al, Nucl. Instr. Meth. A 323 (1992) 162; 
M. Chemarin et al, Nucl. Instr. Meth. A 349 (1994) 345; 
M. Acciarri et al, Nucl. Instr. Meth. A 351 (1994) 300; 
G. Basti et al, Nucl. Instr. Meth. A 374 (1996) 293; 
A. Adam et al, Nucl. Instr. Meth. A 383 (1996) 342 

[8] M. Carena et al, Preprint CERN-TH/99-374 

[9] S.Heinemeyer, W.Hollik, G. Weiglein, Eur. Phys. J. C9 (1999) 343;S. Heinemeyer, W. 
Hollik, G. Weiglein, Phys. Rev. D 58 (1998) 091701;S. Heinemeyer, W. Hollik, G. Weiglein, 
Phys. Lett. B 440 (1998) 296 

[10] M. Carena, M. Quiros and C. E. M. Wagner, Nucl. Phys. B 461 (1996) 407;H. Haber, R. 
Hempfling, A. Hoang, Z. Phys. C 75 (1997) 539 

[11] P. Janot, "The HZHA generator", in "Physics at LEP2", CERN Report 96-01 (1996), 
Version 3, released in December 1999, 
http : // alephwww . cern . ch/~ j anot/Generators . html 

[12] PYTHIA versions 5.722 and 6.1 are used. 

T. Sj6strand, Preprint CERN -TH/71 12/93 (19 93), revised August 1995; Comp. Phys. 
Comm. 82 (1994) 74; Preprint |hep-ph/0001032l (2000) 

[13] KK2f version 4.13 is used; 

S. Jadach, B.F.L. Ward and Z. Was, Comp. Phys. Comm. 130 (2000) 260 

[14] KORALW version 1.33 is used. 

S. Jadach et al, Comp. Phys. Comm. 94 (1996) 216; 
S. Jadach et al, Phys. Lett. B 372 (1996) 289 

[15] KORALZ version 4.02 is used. 

S. Jadach, B.F.L. Ward and Z. Was, Comp. Phys. Comm. 79 (1994) 503 



8 



[16] PHOJET version 1.05 is used. 

R. Engel, Z. Phys. C 66 (1995) 203; 

R. Engel and J. Ranft, Phys. Rev. D 54 (1996) 4244 

[17] F. A. Berends, R. Pittau and R. Kleiss, Comp. Phys. Comm. 85 (1995) 437 

[18] GEANT version 3.15 is used; R. Brun et al, Preprint CERN DD/EE/84-1 (1984), revised 
1987 

[19] H. Fesefeldt, Report RWTH Aachen PITHA 85/02 (1985) 

[20] L3 Collaboration, P. Achard et al, Phys. Lett. B 517 (2001) 319 

[21] S. Bethke et al, Nucl. Phys. B 370 (1992) 310 

[22] A. Read, "Modified Frequentist Analysis of Search Results" in Workshop on Confidence 
Limits, eds. F. James, L.Lyons and Y.Perrin, CERN 2000-05, p. 81 

[23] R. D. Cousins and V. L. Highland, Nucl. Inst. Meth. A 320 (1992) 331 

[24] W. Loinaz and J.D. Wells, Phys. Lett. B 445 (1998) 178;H. Baer and J.D. Wells, Phys. 
Rev. D 57 (1998) 4446;M. Carena, S. Mrenna, C.E.M. Wagner, Phys. Rev. D 60 (1999) 
075010 

[25] OPAL Collaboration, G. Alexander et al, Z. Phys. C 73 (1997) 189. 



9 



The L3 Collaboration: 



P.Achard, O.Adriani? 8 M.Aguilar-Benitez? 5 J.Alcaraz? 5 ' 19 G.Alemanni? 3 J.Allaby? 9 A.Aloisio? 9 M.G.Alviggi? 9 

H. Anderhubf V.P.Andreev?' 34 F.Anselmo? A.Arefiev? 8 T.Azcmoon? T.Aziz? ' 19 P.Bagnaiaf 9 A.Bajof 
G.Baksay? 6 L.Baksay? 6 S.V.Baldew? S.Banerjee? Sw.Banerjee? A.Barczyk? 7 ' 45 R.Barillere? 9 P.Bartalini? 3 
M.Basile? N.Batalova? 4 R.Battiston? 3 A.Bay? 3 F.Becattini? 8 U.Becker? 4 F.Behner? 7 L.Bellucci? 8 R.Berbeco? 
J.Berdugo? 5 PBerges? 4 B.Bertucci? 3 B.L.Betevf M.Biasini? 3 M.Biglietti? 9 A.Biland? 7 J.J.Blaising? S.C.Blyth? 5 
G.J.Bobbink? A.Bohm? L.Boldizsar? 3 B.Borgia? 9 S.Bottai? 8 D.Bourilkov? 7 M.Bourquin? 1 S.Braccini? 1 
J.G.Branson? 1 F.Brochu? J.D.Burger? 4 W.J.Burger? 3 X.D.Cai? 4 M.Capell? 4 G.Cara Romeo? G.Carlino? 9 
A.Cartacci? 8 J.Casaus? 5 F.Cavallari? 9 N.Cavallo? 6 C.Cecchi? 3 M.Cerrada? 5 M.Chamizo? 1 Y.H.Chang? 9 
M.Chemarin? 4 A.Chen? 9 G.Chen? G.M.Chen? H.F.Chen? 2 H.S.Chen? G.Chiefari? 9 L.Cifarelli? F.Cindolo? 

I. Clare? 4 R.Clare? 8 G.Coignet? N.Colino? 5 S.Costantini? 9 B.de la Cruz? 5 S.Cucciarelli? 3 J.A.van Dalen? 1 
R.de Asmundis? 9 P.Deglon? 1 J.Debreczeni? 3 A.Degre? K.Dehmelt? 6 K.Deiters? 5 D.della Volpe? 9 E.Delmeire? 1 
PDenes? 7 F.DeNotaristefani? 9 A.De Salvo? 7 M.Diemoz? 9 M.Dierckxsens? C.Dionisi? 9 M.Dittmar? 7 ' 19 A.Doria? 9 
M.T.Dova? 1 '" D.Duchesneau? B.Echenard? 1 A.Eline? 9 H.E1 Mamouni? 4 A.Engler? 5 F.J.Eppling? 4 A.Ewers? 
PExtermann? 1 M.A.Falagan? 5 S.Falciano? 9 A.Favara? 2 J. Fay? 4 O.Fedin? 4 M.Felcini? 7 T.Ferguson? 5 H.Fesefeldt? 

E. Fiandrini? 3 J.H.Field? 1 F.Filthaut? 1 P.H.Fisher? 4 W.Fisher? 7 I.Fisk? 1 G.Forconi? 4 K.Freudenreich? 7 

C. Furetta? 7 Yu.Galaktionov? 8 ' 14 S.N.Ganguli? P.Garcia- Abia?' 19 M.Gataullin? 2 S.Gentile? 9 S.Giagu? 9 
Z.F.Gong? 2 G.Grenier? 4 O.Grimm? 7 M.W.Gruenewald? 7 M.Guida? R.van Gulik? V.K.Gupta? 7 A.Gurtu? 
L.J.Gutay? 4 D.Haas? R.Sh.Hakobyan? 1 D.Hatzifotiadou? T.Hebbeker? A.Herve? 9 J.Hirschfelder? 5 H.Hofer? 7 
M.Hohlmann? 6 G.Holzner? 7 S.R.Hou? 9 Y.Hu? 1 B.N.Jin? L.W.Jones? Pde Jong? I.Josa-Mutuberria? 5 D.Kafer? 
M.Kaur? 5 M.N.Kienzle-Focacci? 1 J.K.Kim? 3 J.Kirkby? 9 W.Kittel? 1 A.Klimentov? 4,28 A.C.Konig? 1 M.Kopal? 4 
V.Koutsenko? 4 ' 28 M.Kraber? 7 R.W.Kraemer? 5 W.Krenz? A.Kriiger? 6 A.Kunin? 4 P.Ladron de Guevara? 5 
I.Laktineh? 4 G.Landi? 8 M.Lebeau? 9 A.Lebedev? 4 PLebrun? 4 PLecomte? 7 PLecoq? 9 P.Le Coultre? 7 

J.M.Le Goff? 9 R.Leiste? 6 M.Levtchenko? 7 PLevtchenko? 4 C.Li? 2 S.Likhoded? 6 C.H.Lin? 9 W.T.Lin? 9 F.L.Linde? 
L.Lista? 9 Z.A.Liu? W.Lohmann? 6 E.Longo? 9 Y.S.Lu? K.Liibelsmeyer? C.Luci? 9 L.Luminari? 9 W.Lustermann? 7 
W.G.Ma? 2 L.Malgeri? 1 A.Malinin? 8 C.Maiia? 5 D.Mangeol? 1 J.Mans? 7 J.P.Martin? 4 F.Marzano? 9 K.Mazumdar? 
R.R.McNeil? S.Mele? 9,29 L.Merola? 9 M.Meschini? 8 W.J.Metzger? 1 A.Mihul? 2 H.Milcent? 9 G.Mirabelli? 9 J.Mnich? 
G.B.Mohanty? G.S.Muanza? 4 A.J.M.Muijs? B.Musicar? 1 M.Musy? 9 S.Nagy? 6 S.Natale? 1 M.Napolitano? 9 

F. Nessi-Tedaldi? 7 H.Newman? 2 T.Niessen? A.Nisati? 9 H.Nowak? 6 R.Ofierzynski? 7 G.Organtini? 9 C.Palomares? 9 

D. Pandoulas? PPaolucci? 9 R.Paramatti? 9 G.Passaleva? 8 S.Patricelli? 9 T.Paul? 1 M.Pauluzzi? 3 C.Paus? 4 F.Pauss? 7 
M.Pedace? 9 S.Pensotti? 7 D.Perret-Gallix? B.Petersen? 1 D.Piccolo? 9 F.Pierella? M.Pioppi? 3 P.A.Piroue? 7 

E. Pistolesi? 7 V.Plyaskin? 8 M.Pohl? 1 V.Pojidaev? 8 J.Pothier? 9 D.O.Prokofiev? 4 D.Prokofiev? 4 J.Quartieri? 

G. Rahal-Callot? 7 M.A.Rahaman? PRaics? 6 N.Raja? R.Ramelli? 7 P.G.Rancoita? 7 R.Ranieri? 8 A.Raspereza? 6 
PRazis?°D.Ren? 7 M.Rescigno? 9 S.Reucroft? 1 S.Riemann? 6 K. Riles? B.P.Roe? L.Romero? 5 A.Rosca? 

S. Rosier-Lees? S.Roth? C.Rosenbleck? B.Roux? 1 J.A.Rubio? 9 G.Ruggiero? 8 H.Rykaczewski? 7 A.Sakharov? 7 
S.Saremi? S.Sarkar? 9 J.Salicio? 9 E.Sanchez? 5 M.P.Sanders? 1 C.Schafer? 9 V.Schegelsky? 4 S.Schmidt-Kaerst? 
D.Schmitz? H.Schopper? 8 D.J.Schotanus? 1 G.Schwering? C.Sciacca? 9 L.Servoli? 3 S.Shevchenko? 2 N.Shivarov? 2 
V.Shoutko? 4 E.Shumilov? 8 A.Shvorob? 2 T.Siedenburg? D.Son? 3 C.Souga? 4 PSpillantini? 8 M.Steuer? 4 

D. P.Stickland? 7 B.Stoyanov? 2 A.Straessner? 9 K.Sudhakar? G.Sultanov? 2 L.Z.Sun? 2 S.Sushkov? H.Suter? 7 
J.D.Swain? 1 Z.Szillasi? 6 ' 1 X.W.Tang? P.Tarjan? 6 L.Tauscher? L.Taylor? 1 B.Tellili? 4 D.Teyssier? 4 
C.Timmermans? 1 Samuel C.C.Ting? 4 S.M.Ting? 4 S.C.Tonwar? ' 19 J.Toth? 3 C.Tully? 7 K.L.Tung?J.Ulbricht? 7 

E. Valente? 9 R.T.Van de Walle? 1 R.Vasquez? 4 V.Veszpremi? 6 G.Vesztergombi? 3 I.Vetlitsky? 8 D.Vicinanza? 
G.Viertel? 7 S.Villa? 8 M.Vivargent? S.Vlachos? I.Vodopianov? 4 H.Vogel? 5 H.Vogt? 6 I.Vorobiev? 5 ' 28 
A.A.Vorobyov? 4 M.Wadhwa? W.Wallraff? X.L.Wang? 2 Z.M.Wang? 2 M.Weber? PWienemann? H.Wilkcns? 1 
SWynhoff? 7 L.Xia? 2 Z.Z.Xu? 2 J.Yamamoto? B.Z.Yang? 2 C.G.Yang? H.J.Yang? M.Yang? S.C.Yeh? An.Zalite? 4 
Yu.Zalite? 4 Z.P.Zhang? 2 J.Zhao? 2 G.Y.Zhu? R.Y.Zhu? 2 H.L.Zhuang? A.Zichichi?' 19,20 B.Zimmermann? 7 M.Zoller 1 



10 



1 I. Physikalisches Institut, RWTH, D-52056 Aachen, Germany 51 
III. Physikalisches Institut, RWTH, D-52056 Aachen, Germany 5 

2 National Institute for High Energy Physics, NIKHEF, and University of Amsterdam, NL-1009 DB Amsterdam, 
The Netherlands 

3 University of Michigan, Ann Arbor, MI 48109, USA 

4 Laboratoire d'Annecy-le-Vieux de Physique des Particules, LAPP,IN2P3-CNRS, BP 110, F-74941 
Annecy-le-Vieux CEDEX, France 

5 Institute of Physics, University of Basel, CH-4056 Basel, Switzerland 

6 Louisiana State University, Baton Rouge, LA 70803, USA 

7 Institute of High Energy Physics, IHEP, 100039 Beijing, China A 

8 Humboldt University, D-10099 Berlin, Germany 5 

9 University of Bologna and INFN-Sezione di Bologna, 1-40126 Bologna, Italy 

10 Tata Institute of Fundamental Research, Mumbai (Bombay) 400 005, India 

11 Northeastern University, Boston, MA 02115, USA 

12 Institute of Atomic Physics and University of Bucharest, R-76900 Bucharest, Romania 

13 Central Research Institute for Physics of the Hungarian Academy of Sciences, H-1525 Budapest 114, Hungary* 

14 Massachusetts Institute of Technology, Cambridge, MA 02139, USA 

15 Panjab University, Chandigarh 160 014, India. 

16 KLTE-ATOMKI, H-4010 Debrecen, Hungary 11 

17 Department of Experimental Physics, University College Dublin, Belfield, Dublin 4, Ireland 

18 INFN Sezione di Firenze and University of Florence, 1-50125 Florence, Italy 

19 European Laboratory for Particle Physics, CERN, CH-1211 Geneva 23, Switzerland 

20 World Laboratory, FBLJA Project, CH-1211 Geneva 23, Switzerland 

21 University of Geneva, CH-1211 Geneva 4, Switzerland 

22 Chinese University of Science and Technology, USTC, Hcfci, Anhui 230 029, China A 

23 University of Lausanne, CH-1015 Lausanne, Switzerland 

24 Institut de Physique Nucleaire de Lyon, IN2P3-CNRS,Universite Claude Bernard, F-69622 Villeurbanne, France 

25 Centro de Investigaciones Energeticas, Medioambientales y Tecnologicas, CIEMAT, E-28040 Madrid, Spainb 

26 Florida Institute of Technology, Melbourne, FL 32901, USA 

27 INFN-Sezione di Milano, 1-20133 Milan, Italy 

28 Institute of Theoretical and Experimental Physics, ITEP, Moscow, Russia 

29 INFN-Sezione di Napoli and University of Naples, 1-80125 Naples, Italy 

30 Department of Physics, University of Cyprus, Nicosia, Cyprus 

31 University of Nijmegen and NIKHEF, NL-6525 ED Nijmegen, The Netherlands 

32 California Institute of Technology, Pasadena, CA 91125, USA 

33 INFN-Sezione di Perugia and Universita Degli Studi di Perugia, 1-06100 Perugia, Italy 

34 Nuclear Physics Institute, St. Petersburg, Russia 

35 Carnegie Mellon University, Pittsburgh, PA 15213, USA 

36 INFN-Sezione di Napoli and University of Potenza, 1-85100 Potenza, Italy 

37 Princeton University, Princeton, NJ 08544, USA 

38 University of Californa, Riverside, CA 92521, USA 

39 INFN-Sezione di Roma and University of Rome, "La Sapienza", 1-00185 Rome, Italy 

40 University and INFN, Salerno, 1-84100 Salerno, Italy 

41 University of California, San Diego, CA 92093, USA 

42 Bulgarian Academy of Sciences, Central Lab. of Mechatronics and Instrumentation, BU-1113 Sofia, Bulgaria 

43 The Center for High Energy Physics, Kyungpook National University, 702-701 Taegu, Republic of Korea 

44 Purdue University, West Lafayette, IN 47907, USA 

45 Paul Scherrer Institut, PSI, CH-5232 Villigen, Switzerland 

46 DESY, D- 15738 Zeuthen, Germany 

47 Eidgenossische Technische Hochschule, ETH Zurich, CH-8093 Zurich, Switzerland 

48 University of Hamburg, D-22761 Hamburg, Germany 

49 National Central University, Chung-Li, Taiwan, China 

50 Department of Physics, National Tsing Hua University, Taiwan, China 

§ Supported by the German Bundesministerium fur Bildung, Wissenschaft, Forschung und Technologie 

$ Supported by the Hungarian OTKA fund under contract numbers T019181, F023259 and T037350. 

% Also supported by the Hungarian OTKA fund under contract number T026178. 

b Supported also by the Comision Interministerial de Ciencia y Tecnologfa. 

tt Also supported by CONICET and Universidad Nacional de La Plata, CC 67, 1900 La Plata, Argentina. 

A Supported by the National Natural Science Foundation of China. 



11 



Sippd ^ Tin 


MsuSY 


Ma 

iv±2 






Mo 

3 




fGeVl 


[GeVl 


fGeVl 


[GeV] 


[GeVl 


"mh— max" 


1000 


200 


-200 


2M SUSY 


800 


"no mixing" 


1000 


200 


-200 





800 


"large— /z" 


400 


400 


1000 


-300 


200 



Table 1: Values of the MSSM parameters for the three scenarios considered in this Letter. 



"m h — max" scenario 


(m A ,m h ) [GeV] 


(90,90) 


(165,110) 


tan/3 


25 


3 


Preselection 


Data 


2096 


Background 


2044 


e(hA -> bbbb) 


94% 




e(hZ — > bbqq) 


89% 


92% 


hA signal 


7.2 




hZ signal 


0.66 


18.4 


Selection 


Analysis branch 


hA hZ 


hZ 


Data 


25 275 


121 


Background 


28.9 259 


120 


e(hA -> bbbb) 


56% 19% 




e(hZ — > bbqq) 


25% 37% 


63% 


hA signal 


4.3 1.5 




hZ signal 


0.18 0.28 


12.6 


s/b > 0.05 


Data 


18 


42 


Background 


21.3 


42.1 


e(hA -> bbbb) 


56% 




e(hZ — > bbqq) 


28% 


55% 


hA signal 


4.3 




hZ signal 


0.21 


11.0 



Table 2: The number of data, expected background, expected signal events and signal efficien- 
cies in the four-jet channel after preselection, after final selection and after applying a cut on 
the final discriminants corresponding to a signal-to-background ratio greater than 0.05. This 
cut is used to calculate the confidence levels. Numbers are given for two Higgs mass hypotheses 
in the "mh— max" scenario. 
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; mh— max" scenario 



(m A ,m h ) [GeV] 


(90,90) 


(165,110) 


tan/5 


25 


3 


Selection 


Data 


28 


8 


Background 


27.4 


6.4 


e(hA -> bbr+r-) 


41% 




e(hA -> r+r^bb) 


41% 




e(hZ -> bbr+r") 


22% 


33% 


e(hZ — > r+r-qq) 


21% 


32% 


hA signal 


0.50 




hZ signal 


0.03 


0.78 


s/b > 0.05 


Data 


2 


3 


Background 


1.6 


2.8 


e(hA -> bbr+r-) 


34% 




e(hA -> r+r-bb) 


30% 




e(hZ -> bbr+r") 


19% 


30% 


e(hZ — > r+r-qq) 


14% 


29% 


hA signal 


0.39 




hZ signal 


0.02 


0.70 



Table 3: The number of data, expected background, expected signal events and signal effi- 
ciencies in the channels containing tau leptons after selection and after applying a cut on the 
final discriminants corresponding to a signal-to-background ratio greater than 0.05. This cut 
is used to calculate the confidence levels. Numbers are given for two Higgs mass hypotheses in 
the "m h — max" scenario. 
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Expected signal events Expected signal events 



Figure 1: The distribution of data, expected background and expected signal events as a 
function of the logarithm of the signal-to-background ratio, logio(s/b), in the four-jet channel 
for the Higgs boson mass hypotheses a) (mA^h) = (90,90) GeV and b) (mA,mh) = (165,110) 
GeV. Integrated distributions of data and expected background events as a function of the 
expected signal are shown in c) and d). 
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Figure 2: The distribution of data, expected background and expected signal events as a func- 
tion of the logarithm of the signal-to-background ratio, logio(s/b), in the channels containing 
tau leptons for the Higgs boson mass hypotheses a) (mA,mh) = (90,90) GeV and b) {m^^m^) = 
(165,110) GeV. Integrated distributions of data and expected background events as a function 
of the expected signal are shown in c) and d). 
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Figure 3: Distributions of a) logY 56 and b) LAAqq i n the hZ — > AAqq — > qqq / q / q"q" channel. The 
points are data collected at y/s = 203 — 209 GeV, the dashed lines are the expected background 
and the hatched histograms are the signal corresponding to the Higgs boson mass hypothesis 
(mA,m t ) = (30,70) GeV. The signal expectation is calculated within the "no mixing" scenario 
at tan/3 = 0.75 and is multiplied by a factor of 5. 
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Figure 4: Exclusion contours in the (tan/3,mh) and (tan/3, m\) planes at 95% confidence level 
for the "mh— max" and "no mixing" scenarios. The hatched area represents the exclusion 
and the crossed area is not allowed by theory. The horizontally hatched area corresponds to 
niA < 10 GeV and was previously excluded by LEP The dashed line indicates the expected 
exclusion in the absence of a signal. 
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Figure 5: Confidence levels a) (1 — CL mod ) and b) (1 — CL S ) as a function of at tan/5 =15 
obtained for the h boson (solid line), and the H boson (dashed line) searches in the "large— /x" 
scenario. Exclusion contours in the c) (tan/3,mh) and d) (tanP,m^) planes for the "large— /x" 
scenario. The crossed area is theoretically inaccessible, the open area is excluded at 95% 
confidence level and the shaded area is experimentally allowed. The dashed line represents the 
expected boundary of the allowed region in the absence of a signal. 
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